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We study the observational constraints on the Phantom Crossing DGP model. We 
demonstrate that the crossing of the phantom divide does not occur within the framework 
of the original Dvali-Gabadadze-Porrati (DGP) model or the DGP model developed by 
Dvali and Turner. By extending their model in the framework of an extra dimension sce- 
nario, we study a model that realizes crossing of the phantom divide. We investigate the 
cosmological constraints obtained from the recent observational data of Type la Super- 
novae, Cosmic Microwave Background anisotropics, and Baryon Acoustic Oscillations. 
The best fit values of the parameters with lcr (68%) errors for the Phantom Crossing 
DGP model are f2 m , = 0.27l°;°2, P = 0.54±^. We find that the Phantom Crossing 
DGP model is more compatible with the observations than the original DGP model or 
the DGP model developed by Dvali and Turner. Our model can realize late-time accel- 
eration of the universe, similar to that of ACDM model, without dark energy due to the 
effect of DGP gravity. In our model, crossing of the phantom divide occurs at a redshift 
of z ~ 0.2. 

Keywords: Modified gravity; Extra dimensions; Phantom crossing; Cosmic acceleration. 



1. Introduction 

General relativity is the most successful theory of gravity, which passes accurate 
tests in the solar system and laboratories. However, late-time accelerated expan- 
sion of the universe was indicated by measurements of distant Type la supernovae 
(SNIajmS n i s no t possible to account for this phenomenon within the framework 
of general relativity containing only matter and radiation. Therefore, a number of 
models containing "dark energy" have been proposed as the mechanism for the 
acceleration. There are currently many dark energy model s includ ing cosmologi- 
cal constant, scalar field, quintessence, and phantom modeli^^^. However, dark 
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energy, the nature of which remains unknown, has not yet been detected. The cos- 
mological constant is the standard candidate for dark energy. To explain the current 
acceleration of the universe, the cosmological constant must have an incredibly small 
value. But its value cannot be explained by current particle physics, and it is plagued 
with fine-tuning problems and the coincidence problem. 

The mystery of the dark components is based on general relativity, and hence, 
it tells us that what we do not know may be the long distance behavior of gravity 
rather than the energy-momentum components in the universe. In this sense, cos- 
mological observations open up a new window to study the properties of gravity on 
large scales. 

An alternative method for explaining the current accelerated expansion of the 
universe is to extend general relativity to more general theories on cosmological 
scales. Instead of adding an exotic component, such as a cosmological constant to 
the right-hand side (i.e., the energy-momentum tensor) of Einstein's field equation, 
the left-hand side (i.e., the Einstein tensor, which is represented by pure geometry) 
can be modifi ed. Typical models based on this modified gravity a ppro ach are f(R) 
mo del J2QH22] anc [ the Dvali-Gabadadze-Porrati (DGP) modeJ22rEH (f or rev i ewS; 
see Ref.ESJ). 

In f(R) models, the scalar curvature R in the standard Einstein-Hilbert gravita- 
tional Lagrangian is replaced by a general function f(R). By adopting an appropri- 
ate function phenomenologically, f(R) models can account for late-time acceleration 
without postulating dark energy. 

The DGP model is an extra dimension scenario. In this model, the universe is 
considered to be a brane; i.e., a four-dimensional (4D) hypersurface, embedded in a 
five-dimensional (5D) Minkowski bulk. On large scales, the late-time acceleration is 
driven by leakage of gravity from the 4D brane into 5D spacetime. Naturally, there 
is no need to introduce dark energy. On small scales, gravity is bound to the 4D 
brane and general relativity is recovered to a good approximation. 

According to various recent observational data including that of SNld^HiSS it 
is possible that the effective equation of state parameter w e e, which is the ratio 
of the effective pressure p e s to the effective energy density p e g, evolves from being 
greater than —1 (non-phantom phase) to being less than —1 (phantom phase™); 
namely, it has currently crossed —1 (the phantom divide). 

f(R) models that realize the crossing of the phantom divide have been 
studiecp2E3] Qn the other hand, in the original DGP modeP2~125] anc i a p h e _ 
nomenological extension of the DGP model described by the modified Friedmann 
equation proposed by Dvali and TurneJ^, the effective equation of state parameter 
never crosses the w e g = — 1 line. 

In this paper, we explain the "Phantom Crossing DGP model"^^ by further 

1041 

extending the modified Friedmann equation proposed by Dvali and Turner . In our 
model, the effective equation of state parameter of DGP gravity crosses the phantom 
divide line. We investigate the allowed parameter region in detailed comparison with 
recent observational data, and show the validity and the properties of our model. 
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This paper is organized as follows. In the next section, we summarize the orig inal 
DGP model and the modified Friedmann equation proposed by Dvali and TurneJ^l, 
and demonstrate that the effective equation of state does not cross the w e g = — 1 
line in these frameworks. Then, we review "the Phantom Crossing DGP model" 
by additionally extending the modified Friedmann equation proposed by Dvali and 
Turner. We show that the effective equation of state parameter of our model crosses 
the phantom divide line. In Section [3] we study the cosmological constraints on the 



Phantom Crossing DGP model obtained from the recent observations of SNIaPS, 
Cosmic Microwave Background (CMB) anisotropieJ^, and Baryon Acoustic Oscil- 
lations (BAO)^. We show the allowed parameter region and the best fit parameters 
and investigate the validity and the properties of our model. Finally, the conclusions 
are given in Section 01 

2. Model 

2.1. DGP model 

The DGP modeP31 assumes that we live on a 4D brane embedded in a 5D Minkowski 
bulk. Matter is trapped on the 4D brane and only gravity experiences the 5D bulk. 



The action is 



16tt 
1 



M (5) / d b x^/-g (5) R {5) 

bulk 



1fi - M (4) / d 4 x v /=^(i? (4) +L m ), (i) 

1U7T J brane 

where the subscripts (4) and (5) denote quantities on the brane and in the bulk, 
respectively. M(k\ (Mu\) is the 5D (4D) Planck mass, and L m represents the matter 
Lagrangian confined on the brane. The transition from 4D gravity to 5D gravity is 
governed by a crossover scale r c - 

M <4) 

- - J f- (2) 



2^(5) 



On scales larger than r c , gravity appears in 5D. On scales smaller than r c , gravity is 
effectively bound to the brane and 4D Newtonian dynamics is recovered to a good 
approximation. r c is a single parameter in this model. 

Assuming spatial hom ogeneity and isotropy, a Friedmann-like equation is ob- 
tained on the brane^HSl 

H 2 = —p + e-, (3) 

3 r c 

where p is the total cosmic fluid energy density on the brane. e = ±1 represents 
the two branches of the DGP model. The solution with e = +1 is known as the 
self- accelerating branch. In this branch, the expansion of the universe accelerates 
even without dark energy, because the Hubble parameter approaches a constant, 
H = l/r c , at late times. On the other hand, e = — 1 corresponds to the normal 
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branch. This branch cannot undergo acceleration without an additional dark energy 
component. Hence, in what follows, we consider the self-accelerating branch (e = 
+1) only. 

For the second term on the right-hand side of Eq. ([3]) , which represents the effect 
of DGP gravity, the effective energy density is 

3 H (A\ 



and the effective pressure is 



Iff ff 

where H = dH/dt is the differential of the Hubble parameter with respect to the 
cosmological time t. Using Eqs. (@| and ([5]), the effective equation of state parameter 
of DGP gravity is given by 

P r 

w r = — -. (6) 

Assuming that the total cosmic fluid energy density, p, of Eq. ^ contains matter 
and radiation, the effective equation of state of DGP, uv c , can also be exactly 
expressed as 

w rc = , (7) 

1 ~\~ ~ t - 

where il m and f2 r are the normalized energy density of matter and radiation on the 
brane, respectively; i.e., fl m = (8irG/3H 2 )p m and il r — (8irG /3H 2 )p r . (p m oc a~ 3 , 
p r oc a -4 ). 

In realistic ranges of the energy density, il m > and fi r > 0, the value of the 
effective equation of state cannot be less than or equal to —1. That is, the effective 
equation of state never crosses the phantom divide line in the original DGP model. 

2.2. DGP model extended by Dvali and Turner 

1041 

Dvali and Turner^ 1 phenomenologically extended the Friedmann-like equation 
(Eq. ([3])) of the DGP model. This model interpolates between the original DGP 
model and the pure ACDM model with an additional parameter a. The modified 
Friedmann-like equation is ^31 

For a = 1, this agrees with the original DGP Friedmann-like equation, while a = 
leads to an expansion history identical to that of ACDM cosmology. 

Differentiating both sides of Eq. ([8]) with respect to the cosmological time t : we 
obtain the following differential equation: 

2H=-8nG(p + P) + J ^L^, (9) 
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where the dot indicates the derivative with respect to the cosmological time. The 
quantity P is the total cosmic fluid pressure on the brane. 

For the second term on the right-hand side of Eq. ([5J , which represents the effect 
of DGP gravity, the effective energy density is 

3 H a 

and from Eq. ([9]), the effective pressure is 



1 



BttG 



aH 



(r c H) 



(11) 



(12) 



From Eqs. (ITU1) and (jTTJ) , the effective equation of state parameter of the DGP 
model extended by Dvali and Turner is given by 

P a 

W a = . 

Pa 

Fig. Q] shows a plot of the behavior of the effective equation of state of 
the DGP model proposed by Dvali and Turner, w a , versus the redshift z for 
a = 1.00,0.50,0.00,-0.50, and -1.00, assuming fi ro>0 = 0.30, (The subscripts 
designate the present value.) 
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Fig. 1. Effective equation of state, w a , versus redshift z of the DGP model extended by Dvali 
and Turner for a = 1.00,0.50,0.00,-0.50, and -1.00 (top to bottom) assuming Q, mfi = 0.30. For 
a = 1.00, this agrees with the original DGP model, while a = 0.00 corresponds to the ACDM 
model. 
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During the earlier radiation-dominated epoch (z 3> 10 4 ), the effective e quat ion 
of state (Eq. ([T2"])l can also be represented by the following equation with a&*. 

la , , 

w a = -l + —, (13) 

and during the matter-dominated era (10 2 > z^>l), 

«J„ = -1 + |. (14) 

At present, w a is a stationary value close to —1. 

However, when a is positive, the effective equation of state, w a , will exceed — 1 
at all times. For negative a, w a is always less than —1. For a — 0, w a is —1 at all 
times. Based on this analysis, crossing of the phantom divide does not occur in the 
DGP model extended by Dvali and Turner. 



2.3. Phantom Crossing DGP model 

We propose the "Phantom Crossing DGP moder^E] which additionally extends 
the modified Friedmann equation (Eq. ([8])) proposed by Dvali and Turner. Our 
model can realize crossing of the phantom divide line for the effective equation of 
state of the DGP gravity. 

As mentioned in the previous section, the effective equation of state parameter of 
the DGP model proposed by Dvali and Turner, w a , takes the value of greater than 
— 1 for positive a and less than —1 for negative a. When a = 0, w a becomes — 1. 
On the basis of these results, by changing the sign of a (from positive to negative), 
the effective equation of state parameter varies from being greater than —1 to being 
less than — 1. 

Here, we shall assume that the physics that modifies Fricdmann-like equation 
satisfies the following simple requirement: 

• The sign of a in Eq. (j8]) varies from being positive to being negative keeping 
the model as simple as possible. 

In accordance with this requirement, we make the following assumption^! 

a = ft-a, (15) 

where a is the scale factor (normalized so that the present day value is unity) . The 
quantity ft is a constant parameter. In the period when the scale factor a is less 
than the parameter ft (a > 0), the effective equation of state exceeds —1. At the 
point when the scale factor a equals ft, (a = 0), the equation of state's value will 
be —1. In the period when the scale factor a exceeds the parameter ft (a < 0), the 
equation of state will be less than —1. In this way, crossing of the phantom divide 
is realized in our model. 
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Replacing a by (3 — a in Eq. ((5J), the Friedmann-like equation in our model is 
given by 

tt2 8ttG H fj - a 

Differentiating both sides of Eq. ([TBI) with respect to the cosmological time t, the 
following differential equation is obtained: 

2H=-8irG(p + P) + (reg)2 _ (/? _ a) ■ (17) 

For the second term on the right-hand side of Eq. (|16p. representing the effect of 
DGP gravity, the effective energy density is 

3 H 8 ~ a 

Pf> = feG^-(M' (18) 
and from Eq. (|17p. the effective pressure is 

1 



IttG 



(/3 - a)H - aH\n(r c H) , H$- a 



(19) 



Using Eqs. (TTg)) and (TT9"1) . the effective equation of state of our model is given by 

w p = ^. 20 

PB 

We integrate the differential equation (Eq. (fl6|) or Eq. (|T7|) ) in the direction 
from present to past. Hence, the "initial" conditions are given by the parameters 
designating the present value, fi m ,0) (3, and Hq. The value of r c is determined so 
that Eq. (fT6|) is satisfied at the present day. 

Fig. [2] shows a plot of the effective equation of state of our model, wp, versus 
the redshift z (also refer Fig. [3l which shows an enlarged view of this diagram). Our 
model is an extension of the DGP model and realizes crossing of the phantom divide. 
The effective equation of state, wg, of models for (3 = 0.50, 0.25, and 0.10 (assuming 
^m.o = 0.30) crosses the phantom divide line when the redshift z ~ 0.2, 0.8, and 
1.6, respectively. We find that the smaller the parameter (3 is, the older the epoch 
is in which the crossing of the phantom divide occurs. /3 is not necessarily equal to 
the scale factor at the time of crossing the phantom divide, even though Eq. (fl5|) 
is assumed. 

The recent observational data for SNIgP show that crossing of the phantom 
divide line occurs at a redshift z ~ 0.222t221 j n our model, for (3 = 0.50 (when 
^m,o = 0.30), crossing of the phantom divide occurs at z ~ 0.2. 

However, the redshift at the time of phantom crossing is model-dependent. In 
the next section, to test the validity of the Phantom Crossing DGP model, we will 
compare the model and recent observational data in detail. 
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Fig. 2. Effective equation of state of our model, wg, versus redshift z. Red (solid), green (dashed), 
blue (dotted) lines represent the cases of /3 = 0.50,0.25, and 0.10, respectively (assuming f2 mj rj = 
0.30). 
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Fig. 3. Detail of the behavior of wp depicted in Fig. [2] near recent epochs. 
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3. Observational Constraints 

In this section, we study the cosmological constraints on the Phantom Crossing 
DGP model obtained from SNIaPS CMB anisotropics^, and BAcESl observations. 
We carry out detailed investigation of the allowed parameter region, and show the 
validity and the properties of our model. In what follows, we consider the spatially 
flat Universe (k = 0) only. 

3.1. Observational data 
3.1.1. Type la supernovae 

For the SNIa data, we use the nearby+SDSS+ESSENCE+SNLS+HST set of 288 
SNIa by Kessler et al. (2009^311. We adopt the set that used SALT-lfEl for SNIa 
lightcurve-fitting, as the MLCS methocP^ - '^ appears to introduce some systematic 



biases^. The dataset gives the distance modulus at redshift /i h s (zj). 
The theoretical (for a given model) distance modulus is defined by 

fj,(z) = 51og 10 D L +[io, (21) 

where Dl is the Hubble free luminosity distance given by 

*-^£wf>- <22) 

and fj,Q is 

,i = 5 log 10 + 25 = 42.38 - 5 log 10 h, (23) 

h being the Hubble constant Hq in units of 100 km s _1 Mpc . 

We will minimize the statistical x 2 function (which determines likelihood func- 
tion of the parameters) of the model parameters. For the SNIa data, 

288 r i \ i \ i 2 



XSNIa — 



E 

i=i 



(24) 



where is the distance-modulus uncertainty including an additional intrinsic dis- 
persion of crjf* = 0.14 mag^SI. 

Because the nuisance parameter /io (Eq. (|23[1 ) is model-independent, we analyt- 
ically marginalize it as follows: 

b 2 

XsNia = a -— , (25) 



c 

where 

288 



(26) 
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and 



2ns 



-Y — 

We use the XsNia function in combination with that of CMB and BAO data. 



(28) 



3.1.2. Cosmic microwave background 

The CMB shift parameter is one of the least model-dependent parameters extracted 
from the CMB data. Because this parameter involves the large redshift behavior 
(z ~ 1000), it gives a complementary bound to the SNIa data (z < 2). The shift 
parameter R is defined as 



R= y/Q 



m.Q 



Hp 

H(z) 



dz. 



(29) 



where zqmb is the redshift at recombination. Although zcmb depends on the mat- 
ter density VL m q and on the baryon density fij q at the ~ 2 percent level, we fix this 

' ' joy! 

redshift to the 5-year WMAP maximum-likelihood value of zcmb = 1090°'. 

Using the 5-year WMAP data of R obs = 1.710 ± 0.019^2 the \ 2 function for 
CMB is 



XCMB 



R- 1.710 



0.019 



(30) 



3.1.3. Baryon acoustic oscillation 

Observations of large-scale galaxy clustering provide the signatures of the BAO. We 
use the measurement of the BAO peak in the distributio n of luminous red galaxies 
(LRGs) observed in the Sloan Digital Sky Survey (SDSsjSHl i t gives 



0.469 x 



and the parameter A is calculated by 



H 



0.96 
(198 



-0.35 



H(ZBAO) 



1/3 



1 



±0.017, 



zbao Jo 



H{z) 



dz 



2/3 



and zbao = 0.35. The \ 2 function for BAO is 



Xbao 



A - 0.469 x {(0.96 /0.98)- - 35 } 
0.017 



(31) 



(32) 



(33) 



To determine the best value and the allowed region of parameters, we will use 
the maximum likelihood method and need to minimize the following quantity: 



2 _ 2 1 2 1,2 
X - XsNia + Xcmb + Xbao- 



(34) 
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3.2. Numerical results 

We present our main results of constraints from the observational data described in 
the previous subsection. 

In Fig.01 we plot the probability contours in the (f2 m ,o, /3)-plane for the Phantom 
Crossing DGP model. In the red region, crossing of the phantom divide occurs. In 
other words, the effective equation of state of the DGP gravity varies from being 
greater than —1 to being less than —1. The blue and light blue contours show the la 
(68%) and 2tr (95%) confidence limits, respectively, from a combined analysis of the 
SNIa, CMB, and BAO data. We find that much of the la (68%) confidence level 
(C.L.) region is included in the region where the crossing of the phantom divide 
occurs. 
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Fig. 4. Probability contours in the (n m ,o, /3)-plane for the Phantom Crossing DGP model. Cross- 
ing of the phantom divide occurs in the red region. The blue and light blue contours show the la 
(68%) and 2cr (95%) confidence limits, respectively, from a combined analysis of the SNIa, CMB, 
and BAO data. The point denotes the best fit. 



Fig. [5] shows the probability distribution of the energy density parameter of 
matter fl m ^ for the Phantom Crossing DGP model from the combination of SNIa, 
CMB, and BAO data, where parameter /? is marginalized. Fig. [5] shows the proba- 
bility distribution of parameter /3 for the Phantom Crossing DGP model from the 
combination of SNIa, CMB, and BAO data, where parameter Sl m o is marginalized. 
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Fig. 6. ID probability distribution of parameter for the Phantom Crossing DGP model from 
the combination of SNIa, CMB, and BAO data. 
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The best fit values of the parameters with la (68%) errors are 

Q m ,o = 0.27™ (35) 

/3 = 0.54™ (36) 

The best fit value of fl m .o for the Phantom Crossing DGP model is similar to 
the value indicated in the ACDM model. We obtained the comparatively stringent 
constraint for parameter £l m ,o- In Eq. (|15[) . the sign of ft — a varies from being 
positive to being negative by the present day for < ft < 1, because the scale factor 
a is normalized so that the present day value is unity. We find that the 68% C.L. 
region of ft (Eq. (131))) ') is included in the range of < ft < 1. 

In Table[TJ we list the best fit parameters, the y 2 values (Eq. (|34p). and the dif- 
ferences of the Akaike information criterion (AICpSI anc | the Bayesian information 
criterion (BIC)^ , for the ACDM model, the Phantom Crossing DGP model, the 
DGP model proposed by Dvali and Turner, and the original DGP model, from a 
combined analysis of the SNIa, CMB, and BAO data. The definitions of AIC and 
BIC are 

AIC = -21nL + 2fc, (37) 

BIC = -2\nL + k\nN, (38) 

where L is the maximum likelihood, k is the number of free model parameters, and N 
is the number of data points used in the fit. The y 2 value for the Phantom Crossing 
DGP model is the smallest in the four models. We find that the Phantom Crossing 
DGP model is more compatible with the observations than the original DGP model 
or the DGP model developed by Dvali and Turner. However, the values of AIC and 
BIC for the ACDM model are smaller than those for the Phantom Crossing DGP 
model, because the number of parameters of the ACDM model is one less than that 
of the Phantom Crossing DGP model. 



Table 1. Results of observational tests from the combination of SNIa, CMB, and BAO data. 



Model 


Best fit parameters 


x 1 


AAIC 


ABIC 


ACDM 


n m> o = 0.27 


247.262 


0.000 


0.000 


Phantom Crossing DGfESEH 


Q m ,o = 0.27, = 0.54 


246.729 


1.467 


5.137 


DGP by Dvali and TurneiES 


C m ,o = 0.27, a = 0.12 


247.086 


1.824 


5.494 


Original DCF^S 


n ro ,0 = 0.29 


269.002 


21.740 


21.740 



The observational constraints on the DGP model developed by Dvali and Turner 
have been studied in Ref. |2H Although observational data used in our paper are 
not the same as those they used, the result that the best fit value of parameter a 
is closer to than 1 is consistent with that of Ref. US] 

Fig. [7] shows the effective equation of state of DGP gravity (or the cosmological 
constant) versus the redshift z, for the best fit models listed in Table [TJ Only our 
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Phantom Crossing DGP model can realize crossing of the phantom divide line at 
the redshift z ~ 0.2. 

-0.6 
-0.65 
-0.7 
-0.75 
lb, -0.8 
^ -0.85 
-0.9 
-0.95 
-1 
-1.05 

0.2 0.4 0.6 0.8 1 

Z 

Fig. 7. The effective equation of state of DGP gravity (or the cosmological constant), versus the 
redshift z, for the best fit models listed in Table[T] (1) ACDM model, fi TO , = 0.27; (2) Phantom 
Crossing DGP model, O m> o = 0.27, /3 = 0.54; (3) DGP model proposed by Dvali and Turner, 
Q mfi = 0.27, a = 0.12; (4)' Original DGP model, Q mfi = 0.29. 



Fig. [8] shows the distance modulus fi relative to that of a constant expansion 
cosmology /i c , versus the redshift z for the best fit models listed in Table [TJ and the 
combined sample of 288 SNlJ^with the SALT-II lightcurvc fitter. When /i//i c is 
positive, cosmic expansion is accelerating. The Phantom Crossing DGP model can 
realize late-time acceleration of the universe similar to that for the ACDM model, 
without dark energy (and also for the DGP model proposed by Dvali and Turner). 
We will distinguish these by calculating the perturbation in future work. 

Fig.[9]shows the normalized energy density of radiation il r , matter f2 m , and DGP 
gravity Qdgp versus the redshift z in the Phantom Crossing DGP model with the 
best fit parameters in Table [IJ where Hdgp = (87rG/3i7 2 )p,g. pp is the effective 
energy density of DGP gravity defined by Eq. (fTSj) . We find that the universe is 
DGP gravity-dominated near recent epochs. Therefore, in the Phantom Crossing 
DGP model, the late-time acceleration is driven by the effect of DGP gravity not 
by dark energy. 
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Fig. 8. The distance modulus /i relative to that of a constant expansion cosmology fi c , versus the 
rcdshift z for the best fit models listed in Table [l] (the models same as Fig. [7}, and the combined 
sample of 288 SNIa . 




10~ 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 10 5 10 6 

z 



Fig. 9. The normalized energy density of radiation Q r , matter O m , and DGP gravity Qrjcp, 
versus the redshift z in the Phantom Crossing DGP model with the best fit parameter in Table [T] 
(f} m ,0 = 0.27, = 0.54). 
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4. Conclusions 

We demonstrated that crossing of the phantom divide cannot occur in the DGP 
model developed by Dvali and Turner or the original DGP model. We studied the 
Phantom Crossing DGP model by extending the DGP gravity in the framework 
of an extra dimension scenario. Our model realizes crossing of the phantom divide. 
We investigated the cosmological constraints obtained from recent observations. The 
best fit values of the parameters with la (68%) errors for the Phantom Crossing 
DGP model arc n m<0 = 0.27+g;^|, (3 = 0.54±gJ^. We find that the Phantom Cross- 
ing DGP model is more compatible with the recent observations than the original 
DGP model or the DGP model developed by Dvali and Turner. Our model can re- 
alize late-time acceleration of the universe very similar to that of the ACDM model, 
without dark energy due to the effect of DGP gravity. In our best fit model (f2 m ,o = 

0. 27. (3 — 0.54), crossing of the phantom divide occurs at the redshift z ~ 0.2. It is 
known that the self-accelerating branch of solutions in the original DGP model suf- 
fer from ghost-like instabilitie a " '. We will study whether the Phantom Crossing 
DGP model has the ghost problem in future work. 
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